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We investigate the effects of the spin-isospin channel of the Skyrme energy functional on predictions 
for Gamow-Teller distributions and superdeformed rotational bands. We use the generalized Skyrme 
interaction SkO’ to describe even-even ground states and then analyze the effects of time-odd spin- 
isospin couplings, first term by term and then together via linear regression. Some terms affect the 
strength and energy of the Gamow-Teller resonance in finite nuclei without altering the Landau 
parameter g' 0 that to leading order determines spin-isospin properties of nuclear matter. Though 
the existing data are not sufficient to uniquely determine all the spin-isospin couplings, we are able 
to fit them locally. Altering these coupling constants does not change the quality with which the 
Skyrme functional describes rotational bands. 

PACS numbers: 21.30.Fe 21.60.Jz 


I. INTRODUCTION 

Effective interactions for self-consistent nuclear struc¬ 
ture calculations are usually adjusted to reproduce 
ground-state properties in even-even nuclei [I]]. These 
properties depend only on terms in the corresponding 
energy functional that are bilinear in time-reversal-even 
(or “time-even”) densities and currents ||. But the func¬ 
tional also contains an equal number of terms bilinear in 
time-odd densities and currents (see ]|,|| and refs, quoted 
therein), and these terms are seldom independently ad¬ 
justed to experimental data. (For the sake simplicity 
we refer below to terms in the functional as time-even 
or time-odd, even though strictly speaking we mean the 
densities and currents on which they depend.) The time 
odd terms can be important as soon as time-reversal sym¬ 
metry (and with it Kramers degeneracy) is broken in the 
intrinsic frame of the nucleus. Such breaking obviously 
occurs for rotating nuclei, in which the current and spin- 
orbit time-odd channels (linked to time-even channels 
by the gauge symmetry) play an important role. Time- 
odd terms also interfere with pairing, correlations in the 
masses of odd-A and odd-odd nuclei J| y] and contribute 
to single-particle energies and magnetic moments 
|E0| . Finally, the spin-isospin channel of the effective in¬ 
teraction determines distributions of the Gamow-Teller 
(GT) strength. 

The latter are the focus of this paper. We explore the 
effects of time-odd couplings on GT resonance energies 
and strengths, with an eye toward fixing the spin-isospin 
part of the Skyrme interaction. As discussed in our pre¬ 
vious study |ff|, there are many good reasons for looking 


at this channel first. For instance, a better description 
of the GT response should enable more reliable predic¬ 
tions for (3 -decay half-lives of very neutron-rich nuclei. 
Those predictions in turn may help us identify the astro- 
physical site of r-process nucleosynthesis, which produces 
about half of the heavy nuclei with A > 70. 

Our goal is an improved description of GT excita¬ 
tions in a fully self-consistent mean-field model. To this 
end, we treat excited states in the Quasiparticle Random 
Phase Approximation (QRPA), with the residual inter¬ 
action taken from the second derivative of the energy 
functional with respect to the density matrix. This ap¬ 
proach is equivalent to the small-amplitude limit of time- 
dependent Hartree-Fock-Bogoliubov (HFB) theory. We 
proceed by taking the time-odd coupling constants in the 
Skyrme energy functional to be free parameters that we 
can fit to GT distributions. We then check that the cou¬ 
pling constants so deduced do not spoil the description 
of superdeformed (SD) rotational bands. 

Our formulation is nonrelativisitic. In relativistic 
mean-field theory (RMF) the time-odd channels, 

referred to as “nuclear magnetism,” are not independent 
from the time-even ones because they arise from the small 
components of the Dirac wave functions. For rotational 
states, the time-odd effects have been extensively tested 
and shown to be important for reproducing experimen¬ 
tal data (see, e.g., Ref. JIJ]). Only the current terms 
and spin-orbit terms play a role there, however, and the 
time-odd spin and spin-isospin channels of the RMF have 
never been tested against experimental data. 

This paper is structured as follows: In Section || we re¬ 
view properties of the Skyrme energy functional. Section 
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[II reviews existing parameterizations of the functional, 
with particular emphasis on time-odd terms. Our main 
results are in Section IV, where we present calculations 
of GT strength and discuss the role played by the time- 
odd coupling constants. Section ^ describes calculations 
of moments of inertia for selected SD bands. Section VI 
contains our conclusions. We supplement our results with 
six Appendices that provide more detailed information on 
local densities and currents (Appendix ^), early param¬ 
eterizations of time-odd Skyrme functionals (Appendix 

B) , the limit of the infinite nuclear matter (Appendix 

C) , Landau parameters of Skyrme functionals (Appendix 

D) and of the Gogny force (Appendix E|), and the resid¬ 
ual interaction in finite nuclei from Skyrme functionals 
(Appendix |]) . 


challenge is to find an energy functional that incorpo¬ 
rates all relevant physics with as few free parameters as 
possible. The density functional approach as used here 
is equivalent to the local density approximation to the 
nuclear G matrix ||l8|| . 

The energy functional investigated here in detail de¬ 
scribes the particle-hole channel of the effective interac¬ 
tion only. For the treatment of pairing correlations, the 
energy functional has to be complemented by an effec¬ 
tive particle-particle interaction that is constructed in a 
similar way from the pairing density matrix; see pT| for 
details. We use here the simplest functional proportional 
to the square of the local pair density with the coupling 
constants given in [jTT|. 


II. A GENERALIZED SKYRME ENERGY 
FUNCTIONAL 

A. Basics of energy density theory 

Many calculations performed with the Skyrme interac¬ 
tion can be viewed as energy-density theory in the spirit 
of the Hohenberg-Kohn-Sham approach Ea], originally 
introduced for many-electron systems. Nowadays, energy 
density theory is a standard tool in atomic, molecular, 
cluster, and solid-state physics jn|, as well as in nuclear 
physics pi . The starting point is an energy functional £ 
of all local densities and currents p, r, J, s, T, and j that 
can be constructed from the most general single-particle 
density matrix 

P = P{ r, o', r; r', a\ t') = ^2 v l ^(r', T ') r) 

k 

(i) 

(see Appendix |A| for more details), where r, a, and t 
are the spatial, spin, and isospin coordinates of the wave 
function. The Hohenberg-Kohn-Sham approach maps 
the nuclear many-body problem for the “real” highly cor¬ 
related many-body wave function on a system of indepen¬ 
dent particles in so-called Kohn-Sham orbitals ipk- The 
equations of motion for ^ are derived from the varia¬ 
tional principle 

S£ = 0 => hif) k (r,a, t) = e k ip k (r, er,r) , (2) 

where the single-particle Hamiltonian h is the sum of the 
kinetic term t and the self-consistent potential T that is 
calculated from the density matrix 

h = 6 ^ = i+t\p] . (3) 

The existence theorem for the effective energy functional 
makes no statement about its structure. The theoretical 


B. The Skyrme energy functional 

Within the local-density approximation, the energy 
functional is given by the spatial integral of the local 
energy density H(r) 

£ = Jd 3 r H{ r) . (4) 

The energy density is composed of the kinetic term Tfkin, 
the Skyrme energy density 7Yskyrme that describes the 
effective strong interaction between the nucleons, and a 
term arising from the electromagnetic interaction 7f om : 

T~t = kin + Skyrme + T~L em • ( 5 ) 

For the electromagnetic interaction, we take the standard 
Coulomb expression, including the Slater approximation 
for the exchange term. The energy functional discussed 
here contains all possible terms bilinear in local densities 
and currents and up to second order in the derivatives 
that are invariant under reflection, time-reversal, rota¬ 
tion, translation, and isospin rotation [ p 0 | . 

Time-reversal invariance requires the energy density to 
be bilinear in either time-even densities or time-odd den¬ 
sities, so the Skyrme energy density can be separated into 
a “time-even” part 7f even and a “time-odd” part H odd : 

t 

Wskyrme= £ £ + H°? d ) . ( 6 ) 

(=0,1 t 3 =-t 

The sum runs over the isospin t and its third component 
£ 3 . Only the £3 = 0 component of the isovector t = 1 
terms contribute to nuclear ground states and the rota¬ 
tional bands discussed later, while the £3 = ±1 compo¬ 
nents contribute only to charge-exchange (e.g. GT) ex¬ 
citations. In the notation of Refs. fi|.pO| . the time-even 
and time-odd Skyrme energy densities read 

^t™ n = Cf P Ptt 3 + C( Ap Ptt 3 Aptt 3 + O t T Ptt 3 T t t 3 

+ C7 J ptt 3 V • J(t 3 + C/J&,,, 
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( 7 ) 








III. EXISTING PARAMETERIZATIONS 


—Q S S tt 3 + S «3 ' ^ s t*3 + Ct T S tt 3 ' Ttt 3 

+ C t Vs (V • s tt3 ) 2 + Cif tt3 + C? j s „ 3 • V x j tt3 . (8) 

Isospin invariance of the Skyrme interaction makes 
the coupling constants independent of the isospin 21 - 
projection. All coupling constants might be density de¬ 
pendent. Following the standard ansatz for the Skyrme 
interaction, we neglect such a possibility except in C/ 
and Cf. for which we restrict the density dependence to 
the following form 

c p ap o] = cm + (c?[pn m ] - cm) (^)", (9) 

Ct{p o] = C?[0] + (cl [p nm ] - Cl [0]) . (10) 

Here po is the isoscalar scalar density, and p nm is its 
value in saturated infinite nuclear matter. The exponent 
a that specifies the density dependence of C p [po] must 
be about 0.25 for the incompressibility coefficient K 00 to 
be correct |^l] Q. Although this fact does not restrict 
the analogous power in C t s [po], Eq. ©, we keep £ equal 
to a for simplicity here. Usually we will consider energy 
functionals that are invariant under local gauge transfor¬ 
mations which generalize the Galilean invariance of 
the Skyrme interaction discussed in j2) . Gauge invariance 
links three pairs of time-even and time-odd terms in the 
energy functional: 

ci = -ci , C/ = -C t T , C? j =+C? J , ( 11 ) 


These relations fix all orbital time-odd terms, leaving 
only time-odd terms corresponding to the spin-spin in¬ 
teraction with free coupling constants. Relations ( 0 ) 
lead to a simplified form of Eqs. @ -(^), 


t 

H Skyrme = E E [c p t Ptt 3 +ct s 2 tt3 

t=0,l t 3 =-t 


+ 

+ 

+ 

+ 


Ptt 3 Aptt . 3 + C/ As s t t 3 ■ As tt , 3 
Ct ( Ptt 3 T tt 3 — j tt 3 ) U C t (s tt 3 ■ T tt 3 ~ Jtt 3 ) 
C* VJ (Ptt 3 V • J tt 3 + s tt 3 ■ V x j tt3 ) 


c7 s (v-s tt3 ) 5 


( 12 ) 


The time-even terms of the energy functional can be di¬ 
rectly related to nuclear bulk properties such as E/A, the 
saturation density p n . m ., incompressibility, symmetry en¬ 
ergy, surface and surface symmetry energy, and spin-orbit 
splittings. The remaining time-odd terms cannot. 

We will set the coupling constant C t Vs to 0. The term 
it multiplies comes from a local two-body tensor force 
considered in Skyrme’s original papers 0| and discussed 
by Stancu et al. [ p 6 | , but omitted in all modern Skyrme 
parameterizations except the force SL 1 introduced by Liu 
et al. j 2 ?|], which has not been used since. 


The coupling constants of the time-odd Skyrme energy 
functional are usually taken from the (antisymmetrized) 
expectation value of a Skyrme force 0]. When so ob¬ 
tained, the 16 coupling constants of the energy functional 
(0) are uniquely linked to the 10 parameters t i; Xj, Wq, 
and a of the standard Skyrme force (see Appendix || and 
Eq. (p32j)) . Only a few parameterizations rigidly enforce 
these relations, however. Among them are the forces of 
Ref. 0 (e.g. Zfj), SkP |1|, the Skyrme forces of Ton- 
deur |2l], the recent parameterizations SLy5 and SLy7 
and SkX p^]. Most other parameterizations ne¬ 


glect the J 2 term obtained from the two-body Skyrme 
force, setting Cj = 0. Some authors do this for practi¬ 
cal reasons; the J 2 term is time-consuming to calculate, 
and its contribution to the total binding energy is rather 
small. Other authors (see, e.g., |j(J) find that includ¬ 
ing it with a coupling dictated from the HF expectation 
value of the Skyrme force can lead to unphysical solutions 
and/or unreasonable spin-orbit splittings. For spherical 

shapes, the J 2 term contributes to the time-even energy 
density in the same way as the neglected tensor force. 
One might therefore argue that by including the tensor 
force one could counterbalance the unwanted J 2 term 
exactly |3^|. This argument, however, applies neither to 
deformed shapes nor to time-odd fields. Moreover, ne¬ 
glecting this term often violates self-consistency on the 
QRPA level (see below). 

Although one might disagree with the rationale for ne¬ 
glecting the J 2 terms, it is not easy to adjust the coupling 
constants Cf to spectral data. Large values for Cf can 
be ruled out because they spoil the previously obtained 
agreement for single-particle spectra, but there are broad 
regions of values where they influence the usual time-even 
observables too weakly to be uniquely determined [ 0 | . 
Only once in the published literature has there been an 
attempt to do so |3^]. 

All first-generation Skyrme interactions, e.g., SI, SII 
[fidf , and Sill |j(J , used a three-body delta force instead 
of a density-dependent two-body delta-force to obtain 
reasonable nuclear-matter properties. The three-body 
interactions led to a = 1 for Cf in Eq. , but a dif¬ 
ferent density dependence of the C/. a = 1 is too large 
to get the incompressibility Kao right, and causes a spin 
instability in infinite nuclear matter |34| ] and finite nuclei 
fuf (again only within a microscopic potential frame¬ 
work). Both problems are cured with smaller values of 
a (between 1/6 and 1/3 0|) but the second-generation 
interactions that did so still had problems in the time- 
odd channels, giving a poor description of spin and spin- 
isospin excitations and prompting several attempts to de¬ 
scribe finite nuclei with extended Skyrme interactions. 
Krewald et al. |h|], Waroquier et al. [[l7j, and Liu et 
al. |f27| , for example, introduced additional three-body 
momentum-dependent forces. Waroquier et al. added an 
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admixture of the density-dependent two-body delta force 
and a three-body delta force, while Liu et al. considered 
a tensor force. But none of these interactions has been 
used subsequently. 

Van Giai and Sagawa Q developed the more durable 
parameterization SGII, which gave a reasonable descrip¬ 
tion of GT resonance data known at the time and is still 
used today. The fit to ground state properties was made 
without the J 2 terms, however, even though they were 
used in the QRPA. Consequently, in such an approach, 
the QRPA does not correspond to the small-amplitude 
limit of time-dependent HFB. 

All these attempts to improve the description of the 
time-odd channels impose severe restrictions on the cou¬ 
pling by linking them to the HF expectation value of a 
Skyrrne force, leading to one difficulty or another. The 
authors of Refs. [ |18|]3S|l proceed differently, treating the 
Skyrrne energy functional as the result of a local-density 
approximation. The interpretation of the Skyrrne inter¬ 
action as an energy-density functional, besides relaxing 
the restrictions on the time-odd couplings, endows the 
spin-orbit interaction with a more flexible isospin struc¬ 
ture [jfO] [42]] than can be obtained from the standard 
Skyrrne force (43) . Some of the parameterizations used 
here will take advantage of that freedom. But the au¬ 
thors of Ref. Q include only time-odd terms that are 
determined by gauge invariance; the other couplings are 
tentatively set to zero (C t s = C t As = 0). Such a procedure 
is reasonable when describing natural parity excitations 
within the (Q)RPA, but the neglected spin-spin terms are 
crucial for the unnatural parity states that we discuss. 

In this study, we use the energy-functional approach 
( p~2| ) with fully independent time-even and time-odd cou¬ 
pling constants. Our hope is that this more general for¬ 
mulation will improve the description of the GT prop¬ 
erties while leaving the good description of ground-state 
properties in even nuclei untouched. 


IV. GIANT GAMOW-TELLER RESONANCES 

The repulsive interaction between proton particles 
and neutron holes in the J 71 ' = 1 + (spin-isospin) chan¬ 
nel gives rise to a giant charge-exchange resonance in 
all nuclei with excess neutrons. The centroid of the 
resonance (which typically has a width of 5-10 MeV) 
can be roughly parameterized by the simple formula 
E gt -E f = 26A" 1 / 3 - 18.5(iV - Z)A~ 1 , where E F is 
the centroid of the Fermi resonance Q]. This formula, 
however, captures only average behavior; individual cases 
depend on single-particle structure, and in particular the 
spin-orbit splitting. 

The ability to model GT resonances is crucial for pre¬ 
dictions of nuclear /3 decay. Just as the low-lying El 
strength is depleted by the giant dipole resonance, so the 
low-lying GT strength, responsible for /3 decay, is affected 
by the GT resonance. Since one of our future goals is an 


improved calculation of /3-decay rates in nuclei along the 
r-process path, it is important to develop a reliable de¬ 
scription of the GT giant resonance. 


A. Residual interaction in finite nuclei 


Non-self-consistent calculations often use the residual 
Landau-Migdal interaction in the spin-isospin channel: 


5 (rY) 


where Nq is a normalization factor [see Eq. (jD^)[ and k 
and k / are defined in Appendix |^. In most applications, 
only the s-wave interaction with strength g' Q is used, and 
the matrix elements of the force are not antisymmetrized. 
The underlying single-particle spectra are usually taken 
from a parameterized potential, e.g., the Woods-Saxon 
potential. Typical values for g' 0 , obtained from fits to GT- 
resonance systematics, are 1.4 < g' 0 < 1.6 |45|[l?]]. (See 
Ref. [fl8) for an early compilation of data.) Sometimes 
this approach is formulated in terms of the residual in¬ 
teraction between antisymmetrized states. The results 
are similar, e.g., g' 0 = 1.54 in the double-/3-decay calcu¬ 
lations by Engel et al. (ifj. More complicated residual 
interactions, like boson-exchange potentials, have been 
used as well; see, e.g., Refs. [|o| |5^]. Borzov et al. use a 
renormalized one-pion exchange potential in connection 
with a £ = 0 Landau-Migdal interaction of type ([T^) Q . 

A much simpler residual interaction in the GT chan¬ 
nel is a separable (or “schematic”) interaction, v les = 
kqt (<t ■ o - ') (t • t'), where the strength kqt has to 
be a function of A. This interaction is widely used in 
global calculations of nuclear /3-decay . Sarriguren 

et al. |5(]] use it for a description of the GT resonances 
in deformed nuclei with quasiparticle energies obtained 
from self-consistent HF+BCS calculations. They esti¬ 
mate kgt from the Landau parameters of their Skyrrne 
interaction. (The same prescription is used in their cal¬ 
culations of Ml resonances |5^|.) But however useful this 
approach may be from a technical point of view, it is not 
self-consistent. Nor is it equivalent to using the original 
residual Skyrrne interaction; see, e.g., the discussion in 



A truly self-consistent calculation, by contrast, should 
interpret the QRPA as the small-amplitude limit of time- 
dependent HFB theory. The Skyrrne energy functional 
used in the HFB should then determine the residual in¬ 
teraction between unsymmetrized states in the QRPA: 


_ P£_ _ 

Sp(r ll a ll T 1 ;r 2 , cr 2 , t 2 ) 5p{r[, a[, t[; r' 2 , a’ 2 , t£) 

(14) 


The actual form of the residual interaction that con¬ 
tributes to the QRPA matrix elements of 1 + states is 
outlined in Appendix [f| 
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TABLE I. Landau parameters for various Skyrme interac¬ 
tions from relations (B2) and the Gogny forces D1 and Dls. 
Missing entries are zero by construction. 


Force 

9o 

5i 

52 

g'o 

g[ 

52 


SkM* 

0.33 


0.94 


SGII 

0.62 


0.93 


SkP 

-0.23 

-0.18 

0.06 

0.97 

SkI3 

1.89 


0.85 


SkM 

1.77 


0.88 


SLy4 

1.39 


0.90 


SLy5 

1.14 

0.24 

-0.15 

1.05 

SLy6 

1.41 


0.90 


SLy7 

0.94 

0.47 

0.02 

0.88 

SkO 

0.48 


0.98 


SkO’ 

-1.61 

2.16 

0.79 

0.19 

SkX 

-0.63 

0.18 

0.51 

0.53 


D1 

0.47 

0.06 

0.12 

0.60 

0.34 

0.08 

Dls 

0.48 

-0.19 

0.25 

0.62 

0.62 

-0.04 


B. GT strength distributions from existing Skyrme 
interactions 


Before exploring the time-odd degrees of freedom of 
the generalized Skyrme energy functional, we analyze the 
perf ormance of existing parameterizations when relations 
(B2) are used. We examine the forces SkP ]is|| , SGII [fl 8 | . 
SLy4, SLy5 |24j], SkO, and SkO’ |f58| , which all provide 
a good description of ground-state properties but differ 
in details. SkP uses an effective mass m*/?u = 1 and 
is designed to describe both the mean-field and pairing 
effects^. All other forces have smaller effective masses, 
so that m*/m « 0.9 (SkOa;) or even to*/to ~ 0.7 (SGII, 
SLycc). SGII represents an early attempt to get good 
GT response properties from a standard Skyrme force. 
SLy4 and SLy5 are attempts to reproduce properties of 
pure neutron matter together with those of normal nu¬ 
clear ground states. SkO and SkO’ are recent fits that 
include data from exotic nuclei, with particular empha¬ 
sis on isovector trends in neutron-rich Pb isotopes; they 
complement the spin-orbit interaction with an explicit 
isovector degree-of-freedom [Q. All other parameteriza¬ 
tions use the standard prescription C^ J = 3 C^ J ■ 

Residual interactions are often summarized by the 
Landau parameters that appear in Eq. ©• The pa¬ 
rameters can be derived as the corresponding coupling 
constants when Eq. 0 ) is evaluated for infinite spin- 
saturated symmetric nuclear matter (see Appendix O). 
In the literature, the infinite nuclear matter (INM) prop¬ 
erties of the Skyrme interactions are usually calculated 


1 Since the effective mass scales the average density of single¬ 

particle states, it might visibly influence the GT strength 

distribution. 


from Eqs. ( |B2| ). For the generalized energy functional 
(0) discussed here, the time-even INM properties such 
as the saturation density, energy per particle, effective 
mass, incompressibility, symmetry coefficient, and the 
time-even Landau parameters fi, f[ are unchanged, but 
properties of polarized INM and expressions for the time- 
odd Landau parametersand g[ are different. We de¬ 
rive them in Appendix [D|. Here we are most concerned 
with the Landau parameters in the spin and spin-isospin 
channels, 


g 0 = N 0 (2C s 0 +2CZ Pp 1 2 0 /3 ), 

(15a) 

9o = Nq(2Ci + 2Ci /?Pq /3 ), 

(15b) 

9i = -2Aq Cq 0 p 2 J\ 

(15c) 

g'x = - 2 iVo Cl PpT, 

(15d) 


where N 0 is given by ( p5| ) and [3 = ( 37 r 2 / 2 ) 2 / 3 . Val¬ 
ues for some typical Skyrme interactions appear in Ta¬ 
ble |. Higher-order Landau parameters are zero for the 
Skyrme functional 0 )- Some of these values differ from 
those given elsewhere because, unlike other authors, we 
insist on exactly the same effective interaction in the HFB 
and QRPA. The coupling constants Cq and Cl are fixed 
by the gauge invariance of the energy functional, which 
means that Cj — 0 for SGII, SLy4 and SkO, because the 
J 2 term was omitted in the corresponding mean-field fits. 
For these interactions g[ = 0 and g' 0 « 0.9. For SkP and 
SLy5, and SLy7, Cf is relatively large (see Table IV), 
leading to a large g[ ss 1 . 0 , but a cancellation between 
two terms makes g' 0 « 0 . 0 . 

Table | also gives values for the Landau parameters cal¬ 
culated for the Gogny forces D1 [||| and Dls Q from 
the expressions provided in Appendix [^. In the spirit of 
the Gogny force as a two-body potential, one has no free¬ 
dom to choose the time-odd terms independently from 
the time-even ones. (Note that the Gogny force, how¬ 
ever, employs the same local-density approximation for 
the density-dependence as the Skyrme energy functional 
that contributes to the £ = 0 Landau parameters.) The 
higher-order Landau parameters are uniquely fixed by 
the finite-range part of the Gogny force. 

Figures |l] and 0 show the summed GT strength B{GT) 
in 208 Pb and 124 Sn, calculated with all the selected 
Skyrme forces. The ground-state energies are calculated 
as described in Ref. [|ll[, and all strengths are divided 
by 1.26 2 , following common practice, to account for GT 
quenching. Although the GT resonance in 208 Pb comes 
out at about the right energy for SGII, SLy4, SkO, and 
SkO’, it is too low for SkP and SLy5. These latter two 
interactions also leave too much GT strength at small 
excitation energies. It is tempting to interpret these find¬ 
ings in terms of the Landau parameters for these inter¬ 
actions. Schematic models suggest [[Tj that an increase 
of g' 0 results in an increased resonance energy and more 
GT strength in the resonance. The nucleus 208 Pb in¬ 
deed behaves in this way, as can be seen in Fig. l'. The 
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FIG. 1. Summed GT strength in 208 Pb calculated with 
several Skyrme interactions, each corresponding to the Lan¬ 
dau parameters g' 0 and g[ as indicated. The experimental 
resonance energy, taken from Ref. jl(|, is indicated by an ar¬ 



il 5 10 15 20 0 5 10 15 20 

excitation energy (MeV) 

FIG. 2. Same as in Fig. ^ except for 124 Sn. 



-0.2 0.0 0.2 0.4 0.6 0.8 1.0 

go’ 

FIG. 3. Deviation of the calculated GT resonance en¬ 
ergy from experiment, E ca \ c — E ex p t, and fraction of the GT 
strength in the resonance, B rea /Btot, versus Landau parame¬ 
ter g' 0 , calculated for several Skyrme interactions (as indicated 
in the lower right panel) in 90 Zr, 112 Sn, 124 Sn, and 208 Pb. Ex¬ 
perimental values are taken from Ref. Jl6| . 

forces SkP and SLy5, with small values of g' Ql yield more 
low-lying strength and a lower resonance energy than the 
remaining forces which correspond to g' 0 ~ 0.9. 

In 124 Sn, however, this simple picture does not hold, 
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detailed experimental data are from a recent experiment by 
Wakasa et al. 


as Fig. shows. The resonance energies are similar (and 
close to the experimental value) for SkP, SLy5, SkO, and 
SkO’ forces with very different values of g' 0: while SGII 
and SLy4 push the resonance energy too high. Only the 
amount of the low-lying strength seems to scale with g' 0 . 
It is interesting, though, that the related forces SLy4 and 
SLy5 (which predict very similar single-particle spectra, 
but have quite different GT residual interactions) agree 
with the schematic model in that SLy4, with larger g' 0 , 
puts the GT resonance at a higher excitation energy. 

It is clear that the scaling predicted by the schematic 
model is too simple, and Fig. || demonstrates this 
clearly. There we show the calculated strengths B les 
in the GT resonances relative to the sum-rule value 
Btot = 3(IV — Z), and the calculated GT resonance ener¬ 
gies locale relative to the experimental values E expt . [For 
90 Zr, 112 Sn, 124 Sn, and 208 Pb we used E expt = 9.4 MeV, 
8.9MeV, 13.7MeV, and 15.5MeV, respectively p6[ . 
Note that the calculated resonance energy depends on a 
prescription (see [[H]) not strictly dictated by the QRPA.] 
The scatter near g' 0 ss 0.9, in both the resonance energy 
and in the amount of low-lying strength, shows that other 
combinations of parameters in the residual interaction 
besides g' 0 affect the GT distribution. This is not en¬ 
tirely surprising given the compl exity of finite nuclei and 
of the interaction ([hi]). In Sect. IVC we quantify these 
other important combinations and discuss their effects. 

But another factor, this one determined by the time- 


even part of the Skyrme functional, affects the GT dis¬ 
tribution: the underlying single-particle spectrum. Since 
GT transitions are especially sensitive to proton spin- 
orbit splittings, small changes in the time-even part of 
the force can, in principle, move the GT resonance con¬ 
siderably. Sensitivity to the spin-orbit splitting is par¬ 
ticularly obvious in 90 Zr, where detailed information has 
been obtained from a recent experiment by Wakasa et 
al. |6f]]. Unlike in 124 Sn and 208 Pb, which respond to a 
GT excitation in a collective way, the 90 Zr GT spectrum 
is dominated by two single-particle transitions, from the 
neutron lgg /2 state to the proton lgg /2 and I 37/2 states. 
The difference between the locations of the two peaks in 
the GT spectrum is the sum of the proton 1 g spin-orbit 
splitting and a contribution from the residual interaction 
(which can be expected to increase the difference). As 
Fig. 0 shows, all interactions, whatever their value for 
g' 0 , overestimate this difference; the resonance energy is 
always too large, even when the residual interaction is 
switched off completely. 

Most Skyrme interactions give spin-orbit splittings in 



FIG. 5. Relative errors in the spin-orbit splitting (calcu¬ 
lated from the intrinsic single-particle energies) for the forces, 
nuclei, and states indicated. Only splittings between states 
which are both above or both below the Fermi surface are 
included. Other states are affected by core polarization and 
cannot be safely described by the mean held . The forces 
SLy5-SkO are ordered according to their values for g' 0 (see 
Fig. |j). SGII-u and SkM*-u are two recent forces with mod¬ 
ified spin-orbit interactions tailored for future use in GT res¬ 
onance studies |62(]. 
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heavy nuclei that are too large jGl|]. We can therefore 
expect errors in their predicted GT strength distributions 
H0- Figure shows errors in the predicted spin-orbit 
energies for the same forces as in Fig. Interactions such 
as SkI3, SkI4, or SLy4 that overestimate the proton spin- 
orbit splittings give the largest resonance energies (and 
tend to overestimate them). The best interaction, in view 
of the combined information from Figs. and ^], appears 
to be SkO’. Therefore, below, we use its time-even energy 
functional for further exploration of the time-odd terms. 

We have included some new forces in Fig. |j; in a recent 
paper |S2|, Sagawa et al. attempt to improve the spin- 
orbit interaction for the standard Skyrme forces Sill, 
SkM*, and SGII, aiming at better GT-response predic¬ 
tions. They generalize the spin-orbit interaction through 
the condition Cfj 7 ' 7 = — CY J and include the J 2 term 
with a coupling given by Eq. (B2). Although the modi¬ 
fied forces SkM*-u, and SGII-u give slightly better de¬ 
scriptions of GT resonances than the original interac¬ 
tions, they generate unacceptable errors in total bind¬ 
ing energies and do not substantially improve the overall 
description of single-particle spectra in 208 Pb. 

A few remarks are in order before proceeding: (i) The 
spin-orbit splittings shown in Fig. are calculated from 
intrinsic single-particle energies. Since experimental data 
are obtained from binding-energy differences between 
even-even and adjacent odd-mass nuclei, core polariza¬ 
tion induced by the unpaired nucleon, which depends 
partly on time-odd channels of the interaction j^,^], alters 
single-particle energies. The effect is largest in small nu¬ 
clei (of the order of 20% in 16 0), decreasing rapidly with 
mass number [|j. (ii) GT distributions are also affected 
by the particle-particle channel of the effective interac¬ 
tion, but mainly at low energies. The GT resonance is 
not materially altered [jll]], so we can safely neglect the 
particle-particle interaction here. 


C. GT resonances from generalized Skyrme 
functionals 


We turn now to generalized energy functionals in which 
the time-odd coupling constants Cf. C t As , and Cj are 
treated as free parameters; that is, we no longer insist 
that the interaction correspond to a two-body Hamil¬ 
tonian with matrix elements that should be antisym¬ 
metrized. As we showed in Sect. IV B, values of the 


Landau parameter g' 0 alone are insufficient to link the 
properties of the GT resonance to the coupling constants 
of the energy density functional. In this section, using the 
time-even functional of SkO’, we study the dependence 
of the GT resonance on several other combinations of the 
coupling constants as well. Because the isoscalar time- 
odd terms do not affect the GT transitions, we focus here 
on the isovector coupling constants Cf, C As , and Cj. 



FIG. 6. Deviation of calculated and experimental GT 
resonance energies (lower panel) and a fraction of the GT 
strength in the GT resonance (upper panel) for 112 Sn, 124 Sn, 
and 208 Pb, calculated with SkO’ and a modified residual 
spin-isospin interaction. C[ is kept at the Skyrme-force value 
and C As is set to zero. Cf is chosen to be density-independent 
and varied to get 0 < g' 0 < 2.4. g[ = 0.19 in all cases. 


1. Study of C{[pu 


We begin with the simplest case, assuming that (i) the 
functional is gauge invariant, (ii) all time-odd coupling 
constants are density-independent, and (iii) the spin- 
surface term can be neglected, i.e., C As = 0. The only 
remaining free parameter in the spin-isospin channel is 
Cf, whic h is directl y re lated to the Landau parameters 
via Eqs. (|l5b|) and (15d): 


C ‘ “ 2 N 0 


(16) 


where g[ is fixed by Cf [also in Eq. (15d)]. Figure || shows 
results for the GT resonance energy when g' 0 is system¬ 
atically varied from its SkO’ value by altering Cf. We 
have chosen only nuclei that can be expected to exhibit 
a collective response to GT excitations. Non-collective 
contributions may show up, however, when the coupling 
constants are changed. In 124 Sn, for example, a state 
below the resonance collects a lot of strength for small 
values of g' 0 . Only by increasing g' 0 does one push that 
strength into the resonance. Similarly, in 112 Sn a state 
about 5 MeV above the resonance increasingly collects 
strength as g' 0 grows. 

As the underlying single-particle spectra are the same 
for all the cases in Fig. the differences are due entirely 
to the value of g' 0 . With increasing g' Q , the resonance 
energy increases and more strength is pushed into the 
resonance. The increase of E res is nearly linear, but the 
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r (fm) 

FIG. 7. Spatial dependence of C°[p] for various values of 
x=C{ [0]/Cf [pnm], cf. Eq. (fic|) , and g' Q fixed at 1.2. The value 
x = 1 corresponds to no density dependence. For larger values 
of x, the residual interaction becomes more repulsive outside 
the nucleus than inside. When x = 0, Cf\p] vanishes at large 
distances, and for negative values of x, the residual interaction 
becomes attractive outside the nucleus. The density profile 
p(r) used in this plot corresponds to 208 Pb. 

lines for different nuclei have different slopes. It is grat¬ 
ifying that the curves for E ca- i c — C exp t all have a zero 
around the same point, g' 0 ~ 1.2. This value is much 
smaller than the empirical value g' Q « 1.8 derived earlier 
IrilfTTpil for at least two reasons: (i) the influence of 
the single-particle spectrum, and (ii) the inclusion in the 
residual interaction of a p- wave force characterized by g[. 
The latter means that g' 0 = 0 does not correspond to a 
vanishing interaction in the spin-isospin channel. 

2. Study of Cf [0]. 

Thus far we have chosen not to let Cf depend on 
the density. Little is known about the empirical den¬ 
sity dependence of the time-odd energy functional, and 
time-odd Landau parameters calculated from a “realis¬ 
tic” one-boson exchange potential in DBHF show only 
a very weak density dependence Q. Because the ki¬ 
netic spin term C^Si^ • Ti t3 , when evaluated in INM, 
also contributes to the density dependence of the Lan¬ 
dau parameters, the density-dependence of that term 
must either be small or nearly canceled by other time- 
odd terms. In any event, in the following, we investigate 
what happens when Cf depends on the (isoscalar) den¬ 
sity in the “standard” way (|l0|). All nuclei we look at 
have finite neutron excess, which means that the central 
density should be slightly smaller than p nm . 

If g' 0 and g[ are fixed in saturated INM, there is one 
free parameter, C f [0] , with which one can vary the den¬ 
sity dependence (|lC)[). (Cf[p nm ] is fixed by the value 
So[/°nm]=l-2, and we set the exponent £=0.25, as it is 
in the time-even energy functional SkO’.) We continue 
here to assume that gauge invariance holds, and that 

Cf s = 0 . 



strength in the resonance when the ratio x=C{\f)]/Cf [p nm ] of 
parameters defining the density dependence of Cf[p] in Eq. 
( ^~ob is varied. Symbols and scales are as in Fig. |(| 

We vary the parameter Cf[0] between — Cf[p n m] and 
2Cf[p nm ]. Figure [f] shows the spatial dependence of 
Cf[p\ for several values of the ratio x=Cf [0]/Cf [p nm ]. By 
changing Cf[0], one can change both the GT resonance 
energy and the amount of the low-lying strength, even 
with <?o[pnm] kept constant. As Fig. || shows, an increase 
of Cf[0] for a given g’ 0 has almost the same effect as an 
increase of g' 0 for a given Cf [0]. Thus, the INM Landau 
parameters do not tell the whole story in a finite nucleus. 
Figures |?] and || show that the spin-spin coupling has the 
largest effect on the GT resonance when it is located at 
or even slightly outside the nuclear radius. 

3. Study o/Cf s . 

The term Cf is Si t3 • Asi t3 is sensitive to spatial varia¬ 
tions of the isovector spin density. Unlike its (isoscalar) 
time-even counterpart C^ p poApo, it should not be called 
a “surface term” because the spatial distribution of Si is 



FIG. 9. Variation of the GT resonance energy and the 
strength in the resonance when C± s is varied. Symbols and 
scales are as in Fig. || 


9 











































































determined by a few single-particle states that do not 
necessarily vary the most at the nuclear surface. In dis¬ 
cussing the effects of this term, we continue to fix Cf at 
its Skyrme-force value via gauge invariance and choose 
Cf to be density-independent and fixed from Eq. ( 0 ) 
with 3o[Pnm] = 1.2. We then vary Cf s over the range 
of ±30MeVfm°, covering the values obtained from the 
original Skyrme forces. As seen in Fig. |)|, an increase 
of Cf s by 30 MeV fm 5 has nearly the same effect on the 
GT resonance energies as a decrease of g' 0 by 0.2, again 
demonstrating that the value of g' 0 does not completely 
characterize the residual interaction in finite nuclei. A 
new feature of Cf s , apparent from the curves for 112 Sn 
and 208 Pb in Fig. is the ability to move the resonance 
around in energy without changing its strength. 


4 ■ Study of Cf . 

Finally, we investigate the influence on the GT 
strength distribution o f the term Cf S\t 3 • T h 3 , which de¬ 
termines g[ [see Eq. As this term is linked by 

gauge invariance 0 ) to the time-even J 2 term, a fully 
self-consistent variation of Cf would require refitting the 
whole time-even sector of the Skyrme functional. (Note 
that our approach removes the constraints ( |B2| ) that link 
Cf to the time-even coupling constants Cf and C t Ap . 
The constraint was retained, however, when SkO’ was 
constructed.) We leave that task for the future, using 
a gauge-invariance breaking-energy functional here with 
Cf f£_C\ to obtain constraints on Cf for future fits. Fig¬ 
ure |nj shows the change in the GT resonance when g[ is 
varied in the range — 1 < g[ < 1. Increasing gf increases 
the energy of the GT resonance for a given g' 0 . Changing 
g[ by 0.2 has nearly the same effect on the GT resonance 
energy as changing g' 0 by 0.2. (This means that g' Q = 1.2, 


g[ — 0.2, as used here, is consistent with the lower end of 
the values 1.4 < g' 0 < 1.6, g[ = 0.0 given in |I|-|4|].) As 
the curves for 208 Pb and 112 Sn demonstrate, however, the 
amount of strength in the resonance does not necessarily 
change when g[ is varied. 


D. Regression analysis of the GT resonances 


In the previous subsection we explored the dependence 
of the GT resonance energies and strengths on particu¬ 
lar time-odd coupling constants of the Skyrme functional 
while keeping the other coupling constants fixed. These 
results show that the GT properties depend on all the 
coupling constants simultaneously, and the effect of vary¬ 
ing one coupling constant may be either enhanced or can¬ 
celled by a variation of another one. In such a situation, 
linear-regression is needed to quantify the influence of the 
coupling constants. 

We analyze the situation by supposing that the GT en¬ 
ergies and strengths are linear functions of four coupling 
constants, i.e., 

-®reg = e o + eiCi s [0] + e 2 Cf[p nm \ + e 3 Cf s + e^Cf , (17a) 
B G ff =b 0 + hCf[ 0] + b 2 Cf[p nm ] + b 3 Cf s + biC/f. (17b) 


In our linear regression method, the coefficients e ? ; and bi 
are d eterm ined by a least-square fit of expressions (17a) 
and ( [l 7b| ) to the given sample of N calculated QRPA 
results, E G ff c (n) and B^f c (n), n = 1,..., N. The calcu¬ 
lated BfJ (; (n) values have been quenched by the usual 
factor of 1.26 2 . 

The sample of QRPA calculations covers the physically 
interesting range of values for the coupling constants. We 
present here results from a sample defined by the hyper¬ 
cube 



FIG. 10. Variation of the GT resonance energy and the 
strength in the resonance when Cf (and thus g[) is varied. Cf 
is readjusted for each value of Cf so that g' 0 = 1.2. Symbols 
and scales are as in Fig. ^j. 


ffo = 0.6(0.2)1.8, 

7 values, 

(18a) 

Cm/C S 1 [ P nm}= -1(1)2, 

4 values, 

(18b) 

Cf s = - 40(20)40, 

5 values, 

(18c) 

Cf = -40(10)0, 

5 values, 

(18d) 


i.e., for the sample of N = 700. We use g' 0 instead of 
Cf [0] for the regression analysis to avoid combinations of 
the coupling constants t hat lea ve g' 0 too far from 1.2, the 
value advocated in Sect. IVC . 

The left panels in Figs. |11| and 0 contain histograms 
of deviations 


6E GT (n ) = E G J c (n) - E G f , (19a) 

6B GT (n) = B^f c (n) — (19b) 

between the calculated and fitted energies and strengths 
in 112 Sn, 124 Sn, and 208 Pb. The widths of these distribu¬ 
tions illustrate the degree to which the linear regression 
expressions ( 0 ) are able to describe the results of the 
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- 2-10 1 2 - 2-10 1 2 
6E GT (MeV) <® GT (MeV) 

FIG. 11. Distribution of differences between calculated 
GT resonance energies and those from the regression analy¬ 
sis, with all points from the sample (left panels) and with a 
reduced sample (right panels). 
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FIG. 12. Same as in Fig. |ll] except for the strength in the 
GT resonance. 



112 Sn 

124 Sn 

208 Pb 

eo 

5.58100 

10.16000 

8.19600 

ei 

0.00305 

0.00488 

0.00674 

ei 

0.02696 

0.03981 

0.06099 

e3 

0.01690 

0.03474 

0.05897 

e4 

-0.01189 

-0.01767 

-0.02198 

( 7E 

0.0479 

0.199 

0.114 

(JE (full) 

0.585 

0.865 

0.520 

TABLE II. 

Coefficients 

ei obtained by the regression anal- 

ysis, Eq. (17a), of the QRPA GT resonance 

energies in the 

reduced sample (see text). 

The standard deviations o\e from 

this sample are compared to the <r.E(full) from the full sample 
of A=700 points. 


112 Sn 

124 Sn 

208 Pb 

bo 

21.86000 

11.75000 

80.62000 

bi 

0.00963 

0.03107 

0.03079 

62 

0.06143 

0.25120 

0.18620 

b 3 

0.09146 

0.30510 

0.33850 

64 

-0.01730 

-0.10280 

-0.09688 

Ctb 

0.608 

2.10 

2.74 

CTB (full) 

5.20 

3.44 

3.82 


TABLE III. Same as in Table [n] except for the strengths 
of the GT resonances. 


QRPA calculations. One can see that the fit GT reso¬ 
nance energies are generally within about ±1 MeV of the 
calculated ones, and the fit strengths within about ± 6 . 

As can be seen from Figs. |], |s[ |j, and [h], the depen¬ 
dence of both the resonance energy and the strength in 
the resonance on the coupling constants is not linear for 
the entire region of coupling constants. Although our 
sample is restricted to the area around the reasonable 
values, at times we leave the region where the regres¬ 
sion can safely be performed. Furthermore, for certain 
combinations of the coupling constants (especially at a 
weak coupling), there are competing states that carry 
strength similar to that of the GT “resonance.” (These 
states often merge into the resonance at a larger cou¬ 
pling.) Finally, the resonance can be fragmented into 
many (sometimes up to 15) states. Therefore, we remove 
certain areas of parameter space where the determination 
of either the energy or the strength of the GT resonance 
is ambiguous. Such areas are almost always singled out 
by particularly large deviations from the fitted values. 

After reducing the sample in this way, we obtain the 
histograms in the right panels of Figs. [ll] and [1^. These 
illustrate the quality of the regression fits obtained for 
samples of N = 542, 664, and 618 in 112 Sn, 124 Sn, and 
208 Pb, respectively. Tables [H] and HI list the correspond¬ 
ing values of the regression coefficients, as well as the 
standard deviations for the GT energies and strengths 
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within each of the samples. 

Figures [ll| and |lj and the standard deviations ob¬ 
tained in the reduced and full samples (Tables |nj and 
[fill ) show that the description obtained by removing a 
small number of points beyond the region of linearity is 
quite good. The GT resonance energies are now repro¬ 
duced within about ±200 keV or less. The description of 
the resonant GT strengths is also improved, especially in 
112 Sn, although here the linear regression cannot work 
too well because the strengths saturate at strong cou¬ 
pling. Nevertheless, the coefficients listed in Tables [H] and 
III allow a fairly reliable estimate of the QRPA values for 
any combination of the coupling constants. The values of 
the coefficients in Tables |f| and HI show that Cf [p nm l and 


cl 


= —9.172 MeV fm 


( 20 ) 


Ci strongly influence properties of the GT resonance, 
and that both the energies and the resonant strengths 
increase when these coupling constants increase. Cf has 
a weaker effect in the opposite direction, while Cf[ 0 ] is 
less important still. 

Without presenting detailed results, we report here on 
two other attempts at regression analysis. We tried to 
analyze the results for all the three nuclei, 112 Sn, 124 Sn, 
and 208 Pb, simultaneously by adding terms e${N — Z) 
and 65 (N — Z) to the regression formulas ©• Linear 
scaling might be obtained by analyzing the QRPA re¬ 
sults for very many nuclei, where the effects due to shell 
structure could average out. In our small sample, shell 
structure is obviously important. We also tried the re¬ 
gression analysis with CT = —9.172 MeV fm 5 fixed at its 
SkO’ value (see Table |TV|) , and without the terms e 4 and 
&4 in the regression formulae ([l7]), so that the functional’s 
gauge invariance was preserved. The results were not sig¬ 
nificantly different from those when was allowed to 
vary freely. Consequently, our analysis does not allow us 
any constraints on Cg that might be used in future fits 
of the time-even part of the energy functional. 

For the SkO’ coupling 
we obtain U®J = 8.3 MeV, 

112 Sn, 124 Sn, and 208 Pb, 


constants (see Table IV), 


14.2 MeV, and 14.2 MeV in 
and 2 U 8 Pb, and B%g = 100 in 208 Pb. 
These values are close to the corresponding experimen¬ 
tal data: 8.9 MeV, 13.7MeV, 15.5 MeV, from Ref. @, 
and ~ 3(N — Z)/1.26 2 . It is not possible, however, to 
find values of the four coupling constants that reproduce 
these four experimental data points exactly. The rea¬ 
son is that the matrix of corresponding regression co¬ 
efficients is almost singular, resulting in absurdly large 
values of the coupling constants. Clearly a determina¬ 
tion of the coupling constants from experiment would 
require more data. Although charge-exchange measure¬ 
ments have been made on many nuclei, we require spher¬ 
ical even-even nuclei that are not soft against vibrations. 
To fit the relevant coupling constants to data, we would 
need, at the minimum, the ability to treat deformed nu¬ 
clei. Meanwhile, we can make a simple choice of time-odd 
coupling constants from the analysis in Fig. 0. The values 


C?[0] = Cfl Pn 

c As 


I = 120 MeV fm 3 
= 0 


(see Table IV) give g' 0 = 1.2 and g[ = 0.19, which are in 
accord with the data we discuss. But these values by no 
means constitute a fit and are not unique. 


V. A CONSISTENCY CHECK: 

SUPERDEFORMED ROTATIONAL BANDS 

Another phenomenon in which the time-odd part of 
the Skyrme energy density functional plays a role is the 
high-spin rotation of very elongated nuclei. In this sec¬ 
tion we demonstrate that reasonable values for the spin- 
isospin coupling constants found when analyzing the GT 
strength are consistent with the description of superde- 
formed rotational bands. 

When a nucleus rotates rapidly, there appear strong 
current and spin one-body densities along with the usual 
particle densities that characterize stationary (time-even) 
states. The time-odd densities are at the origin of strong 
time-odd mean fields. There are already many self- 
consistent studies of high-spin states available; see, e.g., 
reviews in Refs. The role and significance of 

the time-odd mean-field terms, however, has not been 
carefully studied. Basic features of high-spin states can 
often be well described by models that use phenomeno¬ 
logical mean fields of the Woods-Saxon or Nilsson type, 
where no time-odd terms are explicitly present in the one- 
body potential. (The time-odd densities are, however, 
present there through the time-odd cranking term.) For 
the Gogny interaction [^9|, or within the standard RMF 
models 0. they cannot be independently modified; the 
Gogny interaction is defined as a two-body force (where 
the time-odd terms show up as exchange terms), while all 
time-odd terms appearing in standard RMF models are 
fixed by Lorentz invariance. Within the Skyrme frame¬ 
work, the time-odd terms in superdeformed rotational 
states were analyzed in an exploratory way in Refs. ||^o) . 

Unlike the GT response, rotational bands are influ¬ 
enced by both isoscalar and isovector time-odd channels 
of the effective interaction. In fact, the large effects of 
time-odd coupling constants found in || are mainly due 
to the isoscalar channel; the isovector channel induces 
corrections that are smaller, though non-negligible. The 
SkO’ Skyrme parameterization, which we use for GT cal¬ 
culations, is unstable when the original parameters from 
Eq. (B2) are used in the isoscalar spin channel because 
go < — 1 (a fact that is related to the unusually high value 
of g\ in Table |). This leads to unphysical ferromag¬ 
netic solutions where all spins align when the nucleus 
is cranked. Of course, the value of go does not influ¬ 
ence the GT calculations for even-even nuclei presented 
in our study which focuses on the isovector time-odd cou¬ 
pling constants. Consequently, in the following we em¬ 
ploy a simple spin energy functional using the Skyrme 
force value for Cq , setting Cfr 3 = 0 and neglecting den¬ 
sity dependence. We adopt the value go = 0.4 given in 
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FIG. 13. Dynamical moment of inertia J^ 2 ' > in the superdeformed band of 152 Dy calculated with the SkO’ energy density 
functional and modified time-odd coupling constants. In the upper left panel (corresponding to Fig. ^|) all coupling constants are 
chosen to be density independent, Cf is kept at the Skyrme-force value and C^ 3 = 0. In the upper right panel (corresponding 
to Fig. |s|) the density dependence of C{ is varied keeping g' 0 = 1.2. In the lower left panel (corresponding to Fig. |i)j) C^ s is 
varied, while in the lower right panel (corresponding to Fig. [It]) Cj is varied. See text for the choice of isoscalar time-odd 
couplings. 


p2[ (note that a different definition of the normalization 
factor is used there) to fix Cq. 

We perform the calculations in exactly the same way as 
in Ref. |3f by using the code HFODD (vl.75r) described 
in Ref. [38|. We examine 152 Dy, which is a doubly magic 
superdeformed system. Pairing has a minor influence and 
we neglect it. We focus on the dynamic moment of inertia 

J 2 : 


4h 2 
A Ej 

(from experimental data) or 


= 


d z E 

~dP 


du UJ\ — CU2 


( 21 ) 


( 22 ) 


(in calculations). Figure [Id] shows results of calculations 
when one of the four time-odd isovector coupling con¬ 
stants is varied, while the other ones are kept at the 
values mentioned above. Variations of the coupling con¬ 
stants Cf[0], Cq[p nm ], and Cj have little effect on the 
dynamic moments of inertia in 152 Dy. When C^ s is 
varied, the moments change noticeably, but the general 
trend with frequency is still the same. Thus, altering the 
isovector time-odd couplings does not appear to change 


the quality with which we describe superdeformed rota¬ 
tional bands. Of course, a consistent description of both 
the high-spin data and the GT resonance properties over 
a wide range of nuclei will require a much more detailed 
analysis. 


VI. SUMMARY, CONCLUSIONS, AND 
OUTLOOK 

By exploiting the freedom in the Skyrme energy func¬ 
tional, we have taken significant steps towards a fully 
self-consistent description of nuclear ground states and 
the GT response. Along the way, we debunked the no¬ 
tion that the strength and location of the GT resonance 
in finite nuclei is determined entirely by the Landau pa¬ 
rameter < 7 q . Our analysis also shows this parameter to be 
smaller than previous work indicates. 

There are not enough experimental data for spherical 
even-even nuclei to fix the time-odd isovector coupling 
constants; the ability to do calculations in deformed nu¬ 
clei should help there. We could, however, choose values 
that reproduce the data we do analyze, without spoiling 
our description of high-spin superdeformation. Doing a 
lot better may require improving our time-even energy 


13 





















functionals. GT resonance energies and strengths depend 
significantly on spin-orbit splitting as well as the resid¬ 
ual spin-isospin interaction. Until we are better able to 
reproduce single-particle energies, therefore, a fit of the 
time-odd interaction will be tentative. 

We have not considered isoscalar time-odd interac¬ 
tions. The couplings there will be harder to fix because 
there are fewer data on the response, which is not as col¬ 
lective as in the charge-exchange channel. In addition, 
the isovector time-odd terms will play a role in calcula¬ 
tions of isoscalar observables. Though a lot clearly re¬ 
mains to be done, our work can already be put to good 
use. We will, for example, employ the new values for the 
isovector time-odd coupling constants in future calcula¬ 
tions of beta decay and in the observables that tell us 
about the extent of real time-reversal violation in nuclei. 
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APPENDIX A: LOCAL DENSITIES AND 
CURRENTS 

The complete density matrix p(rat, r'a't') in spin- 
isospin space as defined in (Jlj) can be decomposed into 
the sum of scalar p tt3 (r, r') and vector densities s t t 3 (r, r'), 
where the subscripts denote the isospin quantum num¬ 
bers: 

p(rcrT, rW) 

= i Poo( r , r 0 5 T t' + s 0 o(r, r') • n aa , 6 TT > 

+i +i 

+ <W Plt3(U r ')' r rT' + S lt 3 ( r > r ') •<W r TT' 

t3 = —1 t 3 =—l 

(Al) 

The quantities cr GG ' and t* 3 t , are matrix elements of the 
Pauli matrices in spin and isospin space. In terms of 
these, the local density p , spin density s, kinetic density 
r, kinetic spin density T, current j, and spin-orbit tensor 
J are 


The standard two-body Skyrme force is given by [|2|]33|] 

^Skyrme (G, 1 * 2 ) 

= to (1 + X 0 Pa) <K r l - r 2 ) 

+\ ti (1 + Xi P a ) k' 2 (5(ri - t' 2 ) + S( ri - r 2 ) k 2 
+t 2 (1 + x 2 P a ) k' • J(ri - r 2 ) k 
+ g *3 (1 + x 3 P a ) J(ri - r 2 ) p a (^4^) 

+iWo (<Ti + <r 2 ) • k' x i5(ri - r 2 ) k , (Bl) 

where P a = \{ 1 + <xi ■ <r 2 ) is the spin-exchange op¬ 
erator, k= — L(Vi — V 2 ) acts to the right, and 

k' = ^(V^ — V 2 ) acts to the left. Calculating the 
Hartree-Fock expectation value from this force yields the 
energy functional given in Eq. © with the coupling 
constants: 


Co — |*o + jgh Po 

c i = - + zo) - ^t 3 (§ +^ 3 ) 

c o = - 2 *o(| - x o) - ^ 2 * 3(5 - * 3 ) 
ci = -\to~^t 3P z 
Cq = Te *1 + 2*2 (f + * 2 ) 

Cl = -f *1 (| + ^ 1 ) + 1*2 (| + X 2 J 
Cq = Vj - 2 * 1(2 - x i) + 1 * 2(5 + 


Po 


Po 


X 2 , 


Ci = VJ 


16 1 1 16 
_ _ 9 + 1 1 


H - 

(2 + ^ 2 ) 
Cl AP = M*1 (2 + *1) + 52*2 (| 
C o S = 32 *i(§-si) +25*2(5 


X2 j 

X2 
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TABLE IV. Time-odd coupling constants calculated from Eq. (B2) for the Skyrme interactions as indicated. 


Force 

cg[o] 

(MeVfm 3 ) 

C?[0] 

(MeVfm 3 ) 

Cq [/5nm] 

(MeVfm 3 ) 

C[ [pnm] 

(MeVfm 3 ) 

r T 

'-'O 

(MeVfm 5 ) 

c[ 

(MeVfm 5 ) 

pAs 

(MeVfm 5 ) 

c£ a 

(MeVfm 5 ) 

a 

Skll 

695.860 

239.200 

120.190 

99.573 

0.0 

0.0 

192.660 

62.766 

1/4 

SkI3 

84.486 

220.360 

253.180 

113.940 

0.0 

0.0 

92.235 

22.777 

1/4 

SkI4 

44.038 

231.980 

209.030 

104.120 

0.0 

0.0 

124.590 

37.943 

1/4 

SkO 

373.770 

262.960 

41.421 

84.253 

0.0 

0.0 

70.365 

26.590 

1/4 

SkO’ 

277.910 

262.430 

47.082 

84.154 

-104.090 

-9.172 

42.791 

16.553 

1/4 

SkX 

57.812 

180.660 

-35.639 

81.246 

-7.861 

-23.669 

-4.434 

9.514 

1/2 

son 

271.110 

330.620 

61.048 

91.676 

0.0 

0.0 

15.291 

15.283 

1/6 

SkP 

152.340 

366.460 

-31.328 

78.562 

7.713 

-41.127 

-4.211 

9.757 

1/6 

SkM* 

271.110 

330.620 

31.674 

91.187 

0.0 

0.0 

17.109 

17.109 

1/6 

SLy4 

-207.820 

311.110 

153.210 

99.737 

0.0 

0.0 

47.048 

14.282 

1/6 

SLy5 

-171.360 

310.430 

151.080 

99.133 

-14.659 

-65.058 

45.787 

14.000 

1/6 

SLy6 

-201.460 

309.940 

157.050 

100.280 

0.0 

0.0 

48.822 

14.655 

1/6 

SLy7 

-215.830 

310.100 

158.260 

100.640 

-30.079 

-55.951 

49.680 

14.843 

1/6 


Cf 5 = 


r.VJ 

L '0 

cY J 

°0 


c 


Vs 


64 

= -|Wo 
= —\Wq 
= 0 
= 0, 


(B2) 


nine of which are independent. Although in this ap¬ 
proach rjj = 1, many parameterizations of the Skyrme 
interaction set rjj = 0. That violates the interpretation 
of the Skyrme functional as an expectation value of a real 
two-body interaction and removes the rational e fo r cal¬ 
culating the time-odd coupling constants from (B2). For 
Skyrme interactions with a generalized spin-orbit interac¬ 
tion Q, e.g. for SkI3, SkI4, SkO, or SkO’, the spin-orbit 
coupling constants are given by 


C n VJ = -b 4 - \b\ 


— —I/)' 

— 2^4 


(B3) 


The resulting terms in the energy functional again cannot 
be represented as the HF expectation value of a two- 
body spin-orbit potential (see, e.g., Q), again violating 
the assumptions behi nd th e calculation of the time-odd 
coupling constants in (B2). 

As Eqs. (B2) represent the standard approach to the 


time-odd coupling constants, it is worthwhile to take 
a look at the actual values. Table IV compares them 
for several Skyrme forces. None of these parameteri¬ 
zations was obtained from observables sensitive to the 
time-odd terms in the energy functional. Differences 
among the forces merely reflect various strategies for ad¬ 
justing the time-even coupling constants. Values of the 
density-dependent isoscalar coupling constants Cq, either 
at po = 0 or at po = p nm , are scattered in a wide range. 
This is probably one of the main sources of differences in 
the predictions of the forces for time-odd corrections to 
rotational bands. For the SLyx forces, Cq[ 0] is negative, 


which is unusual; most often this part of the isoscalar 
spin-spin interaction is repulsive at all densities. The 
difference will probably cause visible differences in rota¬ 
tional properties whenever the spin density is large at 
the surface. All the forces agree on the isovector cou¬ 
pling constant Cf, especially at the saturation density, 
i.e., C'l [p nm ] « 100 MeV frn 3 . This simply follows from 
the fact that, assuming Eq. (]b^), Cf is proportional to 
the time-even Cq that is fixed from binding energies and 
radii. 


APPENDIX C: INFINITE NUCLEAR MATTER 

1. Introduction 

Homogeneous infinite nuclear matter (INM) is widely 
used to study and characterize nuclear interactions. 
Some INM properties, such as the saturation density, en¬ 
ergy per particle, and asymmetry coefficient, are coher¬ 
ent, and others, such as the incompressibility Koo and the 
sum-rule enhancement factor, are related to excitations 
and can be used as pseudo-observables to compare with 
predictions of nuclear forces. INM properties are also 
often used to adjust the parameters of effective interac¬ 
tions for self-consistent calculations. These properties at 
large asymmetry are key ingredients for the description 
of neutron stars. (See, e.g., Refs. |2^,^] for a discussion 
on the mean-field level.) 

Most papers deal with spin-saturated INM, in which 
the time-odd channels of the interaction discussed here do 
not contribute. Nothing is known about spin-polarized 
INM, which actually may play some role in neutron stars. 
A stability criterion for this exotic system, derived in 
Ref. j7(j|, was even used to adjust the parameters of the 
SLyz forces in Ref. ||23|j24| ]. We do not consider tensor 
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forces in this work. Their contribution to the properties 
of polarized INM were explored, e.g., in Ref. ||7l| . 

2. Degrees of freedom 

The four basic degrees of freedom of homogeneous INM 
are the isoscalar scalar density po, the isovector scalar 
density p \, the isoscalar vector density Sq , and the isovec¬ 
tor vector density si. They can be expressed through the 
usual neutron and proton, spin-up, and spin-down den¬ 
sities in the following way. 


with = ( 67 r 2 ) 1 / 3 /^ 3 . Here, kp is the “average” 

Fermi momentum of the whole system. The kinetic den¬ 
sity in momentum space for each particle species is given 

by 

T qtr (k) = k 2 pqa(k), (C5) 

and the kinetic density in coordinate space is 

Tqo = W fc F.qa = ijPpl* . ( C6 ) 

where (3 = (37t 2 /2) 2 / 3 . Various kinetic densities in the 
spin-isospin space are given by 


Po — Pn| + Pnl + Ppf + p P i, 

Pi = Pnl + Pnl ~ Ppf ~ Ppl : 
s 0 = Pnl ~ Pnl + Ppl ~~ Ppl: 

Si = Pnl — Pnl - Ppl + Ppl- (Cl) 

Similarly, densities of protons and neutrons with spin up 
and down can be expressed as: 

Pnl = \{P0 + Pi + s 0 + Si) = j(l + I T + la —I - Ira) Po, 

Pnl = 4(^0 V Pi So Si) = 4(1 + It la Ira) PO: 

Ppl — 4(^0 Pi T So Si) — ^(1 I T la ^Tfr) PO: 

Ppl = \{P 0 ~~ Pi ~ So + Si) = j(l — It — I <7 + Ira) Po i 

(C2) 

where I T = pi/po is the relative isospin excess, I a = 
so/ Po is the relative spin excess, and I aT = Si/po is the 
relative spin-isospin excess, with —1 < < +1. 

In symmetric unpolarized INM Ii = 0, while in asym¬ 
metric INM pi 7 ^ 0. Polarized INM has So 7 ^ 0, and 
spin-isospin polarized nuclear matter has Si 7 ^ 0. 


3. Fermi surfaces and kinetic densities 


For INM arbitrary asymmetry, the Fermi energy of 
each particle species is different. Finite spin densities Si 
break the isotropy of INM, creating the possibility that 
the Fermi surface will deform |72| . We are mainly in¬ 
terested in INM with small polarization, so we use the 
approximation that all Fermi surfaces are spherical. 

In the mean-field approximation, p qa ( k), the density of 
particles in momentum space with the isospin projection 
q and spin projection < 7 , is 


Pqafe) 


1 for k < fcF.ijer, 
0 for k > kp^qa- 


(C3) 


In asymmetric polarized INM, the relation between the 
Fermi momenta and the isoscalar scalar density reads 

P ° = 3 tt2 = 67/2 51 5Z k Fga ■ ( C4 ) 

q=p,n a = 


3 Q 5/3 7-1(0) 

t 0 = t„ t + r ni + r pT + r pi = 5 /? p 0 ' ity' 

n = Tnl + Tnl - T pT ~ T pl = § /? 

To = Tnl - T ni + Tpl - T p | = | /? p\ ,Z 

Ti = r nT - Tni - T Pl + T Pl = f P Po^T^ > ( C7 ) 

where Fm\ Fm\ Fm \ and Fm^ are functions of the 
relative excesses: 

F £°> = \ [(1 + I T + la + lor)™ + (l + I T ~ I a - /„)"* 

+ (1 -I T +I a - I aT ) m + (1 -I T ~I a + hr) 1 "] 

(C8) 

This is a straightforward generalization of the corre¬ 
sponding definition for asymmetric unpolarized nuclear 
matter given in |23|| . Similarly one defines 

F^n ] = \ [(! + J r + 1° + 4r) m + (l + I T - 4 - U)"* 

-(1 -I T +I„- I GT ) m - (1 -I T -I„+ /„)”*] , 
f£ ] = i [(1 + I T + I a + I ar ) m - (1 + I T - I a - i aT ) m 

+ (1 -Ir+Ia- I*r) m - (l - I T ~ I a + I^r)™] , 

F^ T) = \ [(1 + It + Ia + I*r) m - (l + I T ~ I* ~ Iar)™ 

(1 I T ~b 1(7 hr) 1 " + (1 -I T -I a + 1^] ■ 

(C9) 

For calculations of INM properties, we also need deriva¬ 
tives of these functions. The first derivatives are given 
by 


dF^ 

dFm ] 

dF^ T) 

— mF^ 

— rnr m _ 1 , 

dir 

dir 

dlrr 

OF i 0) 

dFm ] 

dF^ T) 

- =mF' m ) _ 11 

dir 

dlar 

dl a 

dF^ ] 

dF 1 T) 

dF^ T) 

■ = mF^l 4 , 

dir 

dlrr 

dI T 

dF i 0) 

dFm ] 

dF^ 

= 

dI aT 

dl a 

dl r 


while the second derivatives are 


(CIO) 
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32 F ™ < nP (°) 

o t 2 = ~ l ) F m-2, 

d 2 F$ 


dl? 


= m(m — 1) F^\ 2 , 


(Cll) 


for any i,j = r, a, or. Functions of the order to = 0 and 
to = 1 are rather simple: 


F w — 1 
F 1 (0) = 1 


'= 0 . 
^i (i) =A , 


(C12) 


for any i = Some special values Fm (I T ,J CT , I aT ) 

appearing in limiting cases of INM are 

f' 0) (o,o,o) = i , F«(0,0,0) = 0 
f £\ 1, o, 0) = fW( 0,1,0) = F(°)(0,0,1) = 2 m ~ 1 
Fi T) (1,0,0) = FW(0,l,0)=fW(0,0,l) = 2 m_1 
FW(l,l,l) = FW(l,l,l) = 4 m - 1 . (C13) 


^ = 2 4/3 /?| + (Q + CT + C 0 T + Cf)po 

PO L 

+ (Cq + Cf + Cq + Cf) po 


2/3 

Po 

(C15) 


Expressions ( |B2| ) for an antisymmetrized Skyrme force 
imply that Cq + C p + Cq + C{ = 0, and 


^ = 2 4 / 3 /3| [^ + | *2 (1 + x 2 ) Po 
Po L 


2/3 

Po 


(C16) 


The stability of polarized neutron matter for all den¬ 
sities requires X 2 ~ — 1 & so the SLya: interactions 
take X 2 = — 1 (HJ 2 J]. However, from the energy-density- 
functional point of view, the coupling co nsta nts are in¬ 
dependent, and the second term in Eq. (C15) also con¬ 
tributes to the stability condition. 


5. Pressure, Incompressibility and Asymmetry 
Coefficients 


(i) 

while Fm = 0 if A = 0 and one of the other Ij ’s is equal 
to 1, with the last equal to zero. These functions are 
useful when writing down the equation of state and its 
derivatives. 


4. “Equation of state” of asymmetric polarized 
nuclear matter 


In INM Apu-Jr) = As tt:j (r)(r) = j tta (r) = J« 3 (:r) = 
0. We choose pure neutron and proton states, which leads 
to pi t ± 1 = 0, pi := pi t o, and similarly for all other densi¬ 
ties. We take the z axis as the quantization axis for the 
spin, i.e., St >x = St ,y = 0, St := St )Z , and for the kinetic 
spin density T. As discussed in Refs. [f72|, this breaks the 
isotropy of INM, leading to an axially deformed Fermi 
surface, an effect which we neglect. Adding the kinetic 
term, the total energy per nucleon (i.e. the “equation of 
state”) for the energy functional (j?) and (||) is given by 


n 

Po 


3 h 2 O 2/3 zA 3 ) 
5 2m P Po r 5/3 


+(c p + c{ i 2 + c s Q 1 2 + a r 2 T ) po 

+UC T oF$1 + C1I t F$ 

+Cq I*F$ + ClI.rF^) (3 p* /3 . 


(C14) 


For unpolarized INM one has I a = I aT = 0 which recov¬ 
ers the expression given in Ref. |23|] . 

An interesting special case is polarized neutron matter, 
which is discussed in j70| for the Skyrme interactions. A 
stability criterion derived there from the two-body force 
point of view as outlined in Appendix was used to 
constrain the parameters of the SLya; forces 23, 2J]. In 


this limiting case, one has p n f = po, p n ± = p p \ = p P i = 0, 
which is equivalent to I T = I a = I aT = 1 and leads to 


At the saturation point, all first derivatives of the en¬ 
ergy per nucleon have to vanish and all second derivatives 
have to be positive. The first derivative with respect to 
po is related to the pressure, the second derivative with 
respect to pa is related to the incompressibility, and the 
second derivatives with respect to the A is related to 
the asymmetry coefficients. For symmetric matter, the 
first derivatives with respect to the A vanish because the 
energy per nucleon is an even function of all As. The 
pressure is given by 


P =- 


dE 

dV 


2 dH/po 


which gives 


(C17) 


P = 


2 h- a o/o mVV) 

5 2m P Po r 5/3 

+ (C P + C p I 2 + CqII + C{ll T ) p 2 
+ 3 Po ^ {C P o + C P I 2 + C s Q ll + C s iI 2 t ) 

+(3(c t 0 f^+cif£If 


5/3 P (0) 


\r T t _j_ r T t 

1^0 T VI T r /Q lfj 


- o - 1 5/3 ^ ^1 ± 5/3 
iT P ( CT ) T 

■'O r 5/3 T W 1 ^ 5/3 

The incompressibility is defined as 
18P 


8/3 


K = 


2 d 2 H/p 0 

9p 0 - 


Po ' “ dp 2 


(CIS) 


(C19) 


which, for the Skyrme energy functional ( p!4| ) at the 
saturation point (po = p n . m ., I T = I a = I aT = 0) gives 


K n 


= -I{&-^Po)Pp 2 o / 3 F ^ 3 

,dC p 


+ 2pg' 


dpo 


3 d 2 C p 

Pq- 


d 2 Po~ 


(C20) 
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The asymmetry coefficients are: 


a T 


d(j 


CL(7T 


1 d 2 n/ Po 

2 <5/2 / 7 -=J CT =J CTT = 0 

3 1 ^, + (Co + 3C[) po 


/Jp^ + Cfpo, 


1 d 2 H/p 0 

2 / T =/„=/„=() 

I ^ + (Co + 3C(f) po 


/3po /3 + C 0 s p 0 , 


1 d 2 Tt/p Q 

2 d/2 r 


1 

3 


2m 


+ (Co + 3C^) po 


/9pS /3 + C?po. 


(C21) 


(C22) 


(C23) 


Here, a T is the well-known volume asymmetry coefficient 
of the liquid-drop model, and a a and o CTT are its gen¬ 
eralizations to the spin and spin-isospin channels of the 
interaction. At the saturation point, all asymmetry co¬ 
efficients have to be positive. 


APPENDIX D: LANDAU PARAMETERS FROM 
THE SKYRME ENERGY FUNCTIONAL 

A simple and instructive description of the residual in¬ 
teraction in homogeneous INM is given by the Landau in¬ 
teraction developed in the context of Fermi-liquid theory 
[fio| . Landau parameters corresponding to the Skyrrne 
forces are discussed in Refs. |21,^,^6| j38|j73| |. Starting 
from the full density matrix in (relative) momentum 
space p(k(TT(r'r'), the various densities are defined as 


Poo(k) = YYp^' 7 ^' ( D 1 a) 

a t 

Pits ( k ) = Y Y p{^ TcrT ') T r r'> ( Dlb ) 

a t,t' 

§00 (k) = YYpjkarcr'T) <r aa >, (Die) 

a,a' r 

®it 3 ( k ) = Y Y p(kcrro-V') <7aa' T t 3 t ,, (Did) 

<T,C7 / T,T' 


The kinetic densities are given by r tt3 = p tt3 k 2 , 
T tt:j = s tt3 k 2 . The Landau-Migdal interaction is defined 
as 


F(k 1 a 1 T 1 a l 1 T' 1 :k2<T 2 T2<j' 2 T! l ) 

_ _ S*£_ _ 

5p(ki cri n a\ t[ )Sp{k 2 a 2 r 2 ct 2 r' ) 


= /(ki,k 2 ) + / , (ki,k 2 ) ri • t 2 


+S(ki,k 2 ) cr 1 • cr 2 + g\ kr,k 2 ) (oq ■ ■ r 2 ). (D2) 


The isoscalar-scalar, isovector-scalar, isoscalar-vector, 
and isovector-vector channels of the residual interaction 
are given by 


^ ^ <5poo(ki)<5poo(k2) 

(D3a) 

5 2 £ 


f( kr,k 2 ) = 

opit 3 (ki)dpit 3 (k 2 ) 

(D3b) 

5 2 S 

p( ki ,k 2 )- 5goo(kl )^ 0 o(k 2 ) 

(D3c) 

g (ki,k 2 ) - 

Ssn 3 (ki)dsit 3 (k‘ 2 ) 

(D3d) 


Assuming that only states at the Fermi surface con¬ 
tribute, i.e., |ki| = | k 2 1 = fcp, /, f, 9 , and g' depend 
on the angle 9 between ki and k 2 only, and can be ex¬ 
panded into Legendre polynomials, e.g. 


/(ki,k 2 ) 


1 

ivo 


Yf* p *( 0 )- 


(D4) 


The normalization factor Nq is the level density at the 
Fermi surface 


1 

N~o 


TT 2 h 2 
2m* kp 


150-MeV fm 3 

m* 


(D5) 


A variety of definitions of the normalization factor Nq 
are used in the literature and great care has to be taken 
when comparing values from different groups; see, e.g., 
Ref. Q for a detailed discussion. We use the convention 
defined in |l8). The Landau parameters corresponding 
to the general energy functional (j(ij) are 

/ » =jv “( 2 c » + 4 IS w+ ^'’“ + 2 CS ' a Y - 

fo = Nq{2C{ + 2CI0Poq 3 ), 
go = N 0 (2 Cq + 2Cq (3 Poo 3 )) 
g 0 = N 0 (2C s 1 +2CT ppK?), 
fi = —2A1 0 Cq f3 Pqq 3 , 
f[ =-2N 0 Cl P p 2 0 /3 , 

91 = Ppl' 0 \ 

9i = ~2Ao Ci /3poo 3 - (D6) 


Higher-order Landau parameters vanish for the second- 
order energy functional & but not for finite-range in¬ 
teractions as the Gogny force discussed in the next Ap¬ 
pendix. The Landau parameters provide a stability crite¬ 
rion for symmetric unpolarized INM: It becomes unstable 
for a given interaction when either fe, ge, or g' e is less 
than — (2£ + 1). 
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APPENDIX E: LANDAU PARAMETERS FROM 
THE GOGNY FORCE 


The residual interaction in INM from the Gogny force 


VGogny(ri,r 2 ) 

= H (Wi + ~ MiP a P T )e-^- r ^ 2 /^ 

i= 1,2 

+t 0 (l + *0 Pa) Po ( ri t r2 ) ^( r l - r P 


+iWo (<ti + o 2 ) • k' x <5(ri - r 2 ) k (El) 


(see Appendix [§] for the definition of k, k', P a , and P T ) 
has been discussed in (7^]?^]. Evaluating the expressions 
given in j75| for (k,k',q) = (kF,kp,0), one obtains the 
usual Landau parameters 


fi=J 2 [( 4W * + “ 2Hi ~ 

*= 1,2 

+ (-Wi - 2 Bi + 2H Z + 4 M^f] 

+ 3(o f^o (o: + l)(a + 2) po 

ft=Yi~ ( 2H < + M ^f - ( W i - 250 ^°] 

i= 1,2 

+<^0 (1 + 2®o)Po 

9 l=Y i( 2B ' ~ + (“W + 2H P®f] 

i=l,2 

+<ho 2^0 (1 — 2a; 0 )po 

9t = -Y ( M ^ + W ^) + 3*o Po (E2) 

i=l,2 


where 

= ±tt 3 / 2 /4 iVo Seo 


$l‘> = No e- - iAd) 


4<‘> = !^V?Woe-' 


sinh(z) ( —I— 5 - ) — 


3cosh(^) 


with 2 = \pk%t 2. The normalization factor Ao is again 
given by @). 


APPENDIX F: RESIDUAL INTERACTION IN 
FINITE NUCLEI 

Equation (|lj) gives the most general form of the resid¬ 
ual interaction in finite nuclei. Only a few terms con¬ 
tribute to the 1 + isovector excitations of the even-even 


nuclei we are interested in. First of all, only the isovec¬ 
tor densities contribute. Next, the conditions AJ = 1 
and A-7T = 0 between ground state and excited states im¬ 
ply that the only terms in the energy functional that can 
contribute are quadratic in local tensor or vector parity- 
even densities/currents. As can be seen from Table 2 in 
|po| , all possible contributions are time-odd. One finally 
obtains 


3 (ri,r 2 ) 

5 2 £ 


<5si t (ri)<5si t (r 2 ) 

2C® [poo] 5(ri - r 2 ) 

-.A. 


(<Ti ■ <T 2 ) (n • r 2 ) 


+ 3 (Cl - 4C^ s ) (k' 2 <5(n - r 2 ) + <y(n - r 2 )k 2 ) 


+ (3 C[ + 4Cf s ) k' • 6{n - r 2 ) k 


O 1 02 T 1 • T 2 


-2iCy J r i • f 2 (<xi + o 2 ) • k' x 5(ri - r 2 ) k 

where k and k' are defined in Appendix [b|. Since the 
coupling constants depend only on the scalar isoscalar 
density pooj there are no rearrangement terms in the 
spin-isospin channel of the residual interaction. Unsym¬ 
metrized proton-neutron matrix elements of this inter¬ 
action are to be inserted into the QRPA equations as 
outlined in Ref. Ol. 
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